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The temperature-dependent uniform magnetic susceptibility of interacting electrons in one di- 
mension is calculated using several methods. At low temperature, the renormalization group reaveals 
that the Luttinger liquid spin susceptibility \ (T) approaches zero temperature with an infinite slope 
in striking contrast with the Fermi liquid result and with the behavior of the compressibility in the 
absence of umklapp scattering. This effect comes from the leading marginally irrelevant operator, 
in analogy with the Heisenberg spin 1/2 antiferromagnetic chain. Comparisons with Monte Carlo 
simulations at higher temperature reveal that non-logarithmic terms are important in that regime. 
These contributions are evaluated from an effective interaction that includes the same set of di- 
agrams as those that give the leading logarithmic terms in the renormalization group approach. 
Comments on the third law of thermodynamics as well as reasons for the failure of approaches that 
work in higher dimensions are given. 



I. INTRODUCTION 

There are a number of organic conductors, including the Bechgaard salts for example, for Tftt4 cn the distinctive 
behavior of one-dimensional interacting electrons is observed over a wide range of temperature lira Among the char- 
acteristics observed is the asymptotic low-frequency, long-wavelength electronic behavior which, in one dimension, 
belongs to the universality class of Luttinger liquids. This universality class plays in one-dimension the role of Landau 
Fermi liquid theory in higher dimension, providing a framework to understand the .occurence of power law behavior, 
spin-charge separation and various other characteristics of one-dimensional systemsa. 

While much is known about the predictions of Luttinger liquid theory, non-singular quantities, such as the uniform 
magnetic spin susceptibility \, are not completely understood theoretically. Experimentally, \ (T) can be accurately 
measured as a function of temperature using a number of techniques, including Knight shift in Nuclear Magnetic 
Resonance (NMR) experiments. In higher dimension, the theoretical situation for x (T) is clear within Landau Fermi 
liquid theory. The prediction for this key quantity is a Pauli-like susceptibility whose absolute value is enhanced biz. 
interactions. For -Luttinger liquid theory in one dimension pioneering work was done by Dzyaloshinskii and Larkin.cl 
and by Lee et ala A few years ago, BourbonnaisO, using the renormalization group, readdressed this problem. The 
conclusion that emerges from these works is that the temperature dependence of x (T) remains important up to T = 0, 
contrary to the Fermi liquid result in high dimension where temperature dependence shows up only on the scale of 
the Fermi energy. The renormalization group approach is expected to give the correct asymptotic low-temperature 
behavior of the Luttinger liquid. To obtain quantitative confirmation of these results from a specific microscopic 
model, it is customary to study the Hubbard model that has an exact Bethe ansatz solution in one dimension and 
belongs to the Luttinger- liquid .universality class. Despite the exact solution, the Bethe ansatz magnetic susceptibility 
must be computed numericallyfl Additional results have been obtained recently by Jiittner et al.B through a quantum 
transfer matrix method, by Mila and Penc through world line quantum Monte Carlo simulations Jj and by MoukouriliEl 
through Density Matrix Renormalization Group techniques. It is important to notice that the low temperature limit 
is difficult to reach numerically .using .any of the available numerical approaches, including Bethe ansatz. 

In this paper, following Rcfs.El andU, we derive the renormalization group (RG) prediction for the temperature- 
dependent uniform magnetic susceptibility of the g-ology Hamiltonian. We show that the uniform magnetic suscepti- 
bility x {T) of this Luttinger liquid .aoproaches zero temperature with an infinite slope, in analogy with the spin 1/2 
antiferromagnetic Heisenberg chainE3£j. This effect comes from backscattering of right- and left-moving electrons, 
which is the leading marginally irrelevant operator. We show that this phenomenon appears in a temperature range 
where no numerical calculation has been done yet. We also obtain predictions at higher temperature by using nu- 
merical simulations done on the Hubbard Hamiltonian at quarter filling. To obtain quantitative agreement with the 
simulations, we find it necessary to include non-logarithmic contributions that are beyond the RG treatment. This is 
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done by using an approach inspired from the KanamoriE-a-BruecknertL^I theory valid in higher dimensioril3. The subset 
of diagrams that must be resummed in one dimension is suggested by the RG and it does differ from the subset used 
in higher dimension. The reasons for the failure of higher-dimensional approaches are also discussed. 

The </-ology and Hubbard Hamiltonians are introduced in Sec. II. In Sec. Ill, we present the results of numerical 
calculations and in Sec. IV the RG calculation, including the prediction of the infinite slope as T — + 0. Sec. V 
discusses the comparison between numerical results, RG prediction and the reasons for the failure of higher-dimensional 
approaches. We conlude with a summary of our main results and general comments on the range of applicability of 
Luttinger-liquid theory. 



II. HUBBARD MODEL AND CONNECTION WITH G-OLOGY 



The simulations are done for the Hubbard model 

H = -t (4* c j<? + CjaCiv) + U J2 n *T n a C 1 ) 



<l]><7 



where units of energy are chosen such that the hopping matrix element t equals unity in the simulations. The creation 
(annihilation) operators c\ a (ci a ) create (annihilate) electrons of spin a in the orbital located on site i with position r^. 
Only nearest-neighbor hopping is allowed. The last term, with the usual occupation number operators 7ijj = c irj Ci ai 
represents the short-tange repulsion U, felt by the electrons when they occupy the same orbital at site i 



As is well knownjL3 there are logarithmic divergences in the perturbative treatment of t 



Hubbard model in 
L7 ' or most easily by a 



one dimension. These can be handl6d|-dipGctly by infinite resummations of parquet diagrams, 

renormalization group treatmentJaJilM In this case, the Hubbard Hamiltonian is not a fixed-point Hamiltonian. 
It is necessary to consider the renormalization group flows in a more general space of Hamiltonians called g-ology 
Hamiltonian. In the rest of this section, we recall how to cast the Hubbard Hamiltonian as a special case of g-ology. 
First, it is useful to rewrite it in the form 

U 



E ( c L c i ff + 4*°*°) + ~2 E c \a4a> c i<y' c i<y ( 2 ) 



H = —i 

<ij>(T iaa 



which, compared with Eq. (|T|) contains an additional term that can be absorbed in a chemical potential shift. The 
g-ology Hamiltonian being defined in Fourier space, we take 

c a {k) = ±=Y j er^ Cia ; c t (fe) = _L jre ikri 4* (3) 

where we have chosen unity for the lattice spacing so that the number of sites and the system size are both equal to 
L. Using these variables, and neglecting umklapp processes, we can write, 

= ]T(-2icosfc)4 (fc) Cc7 (fc) + I- ]T ci(k)cl(k')c a ,(k' + q )c a (k-q). (4) 

k,a k,k' t q,a,a' 

In recent versions of the renormalization grou the full cosine dispersion relation can be taken into account, but 
in the more usual version that we consider here, the dispersion relation is linearized around the two Fermi points ±fcp, 
and one considers only scatterings around and between these points. To rewrite the Hubbard Hamiltonian Eq.(|^) 
in a way that highlights the processes that are allowed by the Pauli principle, the sum over momentum transfers q 
is divided also in three pieces: q f=a 0, and q w ±2kp. Furthermore, one introduces a lower index p to the creation- 
annihilation operators that, for the moment, just indicates if the allowed particle momenta are mostly around the 
+kp Fermi point (right-moving (+)), or around the — kp Fermi point or (left-moving (— )). This rearrangement gives, 
after one allows the k sums to run from —ko+pkp to ko + pkp with kg a cut-off wave vector of the order of kp . 

H e P ( fc ) C )U ( fc ) c v,<y ( k ) 

k,cr.p 

i ko +pk F k -\-p' k F 

+ E E E E Y. C lAk)cl,, a ,{k')c p ,^{k' + q)c p ,„{k-q) 

0-,f'=T (p,p'=±) (k=-k +pk F ) (k' = -k +p'k F ) <? 
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+ §Z E E E E c +-( fc ) c -^( fc ')c + ,.'(fc' + 9 + 2Mc-, ff (fc- ( z-2fc F ) 

t7,cr'=T (fe=— fco+fcy) (fc' = — fco — fej?) 9 
1 feo+feF k -k F 

+ U E E E E c -,-( fc ) c U'( fc ') c -^( fc ' + '?- 2 Mc + ^(fc-g + 2fc F ) (5) 

<t,ct'=T (fc'=-fc +fe F ) (fe=— fc — fc F ) g 

where the linearized dispersion relation is 

e p (k) = pvf (k — pkp) ; vp = 2isinfcp. (6) 

The restrictions on momentum transfer q should normally be set to avoid double counting various scattering processes. 
However, it is simpler to introduce additional states that linearly extrapolate the right and left-moving electron 
dispersion relations. In otherwords, strikctly speaking we should have c +iCT (k) — c a (k) for < k < tt and c_ !CT (A;) = 
Co- (k) for — 7r < k < while, instead, we add states in such a way k runs from — oo to +oo for both cases p — ±1. Each 
of these sets, p = ±1, arc defined as a "branch" with the corresponding dispersion relation e p (k) = pvf {k — pkp) . 
The added unphysical states should not contribute appreciably because of the large energy denominators, the Pauli 
principle and the cut-offs in the sums over k and k' that regularize perturbation theory. Hence, in the weak to 
intermediate coupling regime only scatterings near the Fermi surface are important so we can assume that the sum 
over q that is left for each of the three pieces is free to run from — oo to +oo. In the g-ology notation, the last two 
terms of Eq.(||) are regrouped into 2k F scatterings with an interaction constant gx, while the first interaction term is 
divided into q ~ scatterings on the same branch (54), and q w scatterings between two different branches (52), 
namely 

H a = E e P W c l,<y c p>° ( fc ) 

i ko+pkp k ~pk F 

+ |E E E E E 4,Ak)cl p ^(k')c p ^(k' + q + 2pk F )c- p , <T (k-q-2pk F ) 

q o-,er' = T P=± (k=-k +pk F ) (k' = -k -pk F ) 
I k +pk F kg— pk F 

+|EE E E E 4,Ak)cl p , a ,W)c- P ,cW + q)c p , t ,(k-q) 

q <t,ct' = T P=± (k=-k a +pk F ) {k' = -k -pk F ) 
I kg+pkp k a +pk F 

+ |ZE E E E E <Ak)cl a ,(k')c p ^(k' + q)c p , a (k-q) (7) 

q ct,ct' = T P=± (k=-k a +pk F ) (k< = -k +pk F ) 

In this notation then, the space of parameters is closed under the renormalization-group (RG) induced flow. 

The Hubbard Hamiltonian and the above g-ology Hamiltonian clearly differ since the dispersion relation is linearized 
and (7-ology contains fewer scattering terms than the original Hubbard model. If one assumes that the high-energy 
processes that were dropped do not influence the Physics, then g\ = g2 = 94 = U for the Hubbard model. It is 
important to notice, however, that this identification is an approximation. In fact the g-ology Hamiltonian H g is a 
small-cut-off limit of the Hubbard model and is strictly related to the Hubbard model only in a RG sense. Generally, 
the appropriate initial values of the coupling constants entering the RG would be such that, for example, 174 ^ U . Here 
by rotational invariance parallel and perpendicular components of the coupling constant g\ are identical, a property 
that is preserved by the RG transformation on the g-ology Hamiltonian. 



III. QUANTUM MONTE CARLO RESULTS 



We are interested in the quarter-filled case, which corresponds to a large class of organic conductors. Earlier 
results for the uniform magnetic susceptibility _x include the following. First, zero temperature results that have 
been obtained from the Bethe ansatz by ShibatJ and are shown for U/t = 2, and U/t = 4 as the left-most points 
on Fig 1. Recently, Jiittner et aln used a new approach based on the Trotter-Suzuki mapping and a subsequent 
investigation of the quantum transfer matrix to obtain the temperature-dependent results shown by the solid line in 
Fig.l. Earlier results had been obtained using a world line algorithm by Mila and Penc.El The corresponding results, 
at U/t — 4 are linked by a dashed line in Fig.l. These results disagree substantially from those of Jiittner et al.U 
As we show below, our Monte Carlo simulations confirm the latter results. The calculations of Ref.u are probably 
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Fig. 1 
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FIG. 1. Monte Carlo simulation results for the temperature dependent susceptibility as defined by Eq.(g|) and T\ = S (q = 0). 
Our Monte Carlo results (2 x 10 5 measurements) are shown by symbols with error bars for different values of U. For [7 = 4 the 
extrapolation of our results to At = is also showipifor three temperatures. They are, within error bars, in agreement wit 
the quantum transfer matrix results of Juttner et al!u shown by the seJid line. Points joined by the dashed line are from Ref. 
while the zero-temperature results shown by symbols are from ShibaE3. 
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less accurate because they were obtained not only from the usual extrapolations to At = and L — oo, they also 
required extrapolation to q = of the finite wave vector results because the world-line algorithm is for the canonical 
ensemble. Finite-temperature DMRG calculations^ cannot yet reach sizes large enough to be compared with the 
results of Fig.l, although for larger values of the interaction it would be possible. .— . 

We have used the so-called determinantal Monte Carlo method (BSS algorithm)c3 to obtain the other finite- 
temperature results in Fig.l for U/t = 2, 3, 4. Since equal-time quantities give better statistics, the quantity that was 
computed is the magnetic structure factor for an L site lattice 

S = \ cxp ( ri ~ r J')] ((™ 1 >T ~ n ^ ( n J'T _ n i.l)> • ( 8 ) 

This quantity at q = is trivially related to the magnetic susceptibility through the fluctuation-dissipation theorem 
(S (0) = T\) ■ The results shown by isolated points with error bars in Fig.l are for a 30 site chain at quarter-filling. 
The units are chosen so that t = 1. About 2 x 10 5 measurements were taken for each point. For U = 4 and three 
temperatures, namely T — 1/3, 1/6 and 1/9, we have done the extrapolation to At = using three values of At. 
These, results are plotted as crosses in Fig.l. Our extrapolated results are in excellent agreement with those of Jiittner 
ct al.Ej, while our unextrapolated results systematically underestimate those of Jiittner et aln by at most a few percent 
at the highest temperatures. We conclude that the discretization step that we used in imaginary time, At = 1/8, 
leads to a systematic underestimation (of ordercj (At) 2 ) of the susceptibility which is smaller than the statistical 
error on the figure except perhaps for the two highest temperatures where the extrapolated results are within two 
error bars of the unextrapolated ones. 

As we shall see in the RG treatment, interactions cause a singularity in the temperature derivative of the q — 
susceptibility at T = 0. Observation of this singularity would require huge system sizes. However, at higher 
temperature, where the susceptibility is regular, finite-size effects should be negligible when the thermal de Broglie 
wavelength £ t h = Vp/ i^T) is smaller than the system size. For system sizes 30 and quarter- filling, this criterion means 
that finite-size effects should be small for T > 0.02. In fact, for the lowest temperature we have considered T « 0.1, 
we have £th ~ 5. We have verified for T > 0.1 that indeed our results are the same, within statistical accuracy, for 
sizes 10,20,30, 40. 

The results shown in Fig.l do not contain the factor 1/2 for each external spin vertex, hence they are larger by a 
factor of four than those defined in the following section. 



IV. RENORMALIZATION GROUP APPROACH 



The renormalization group provides a useful tool to understand the g-ology Hamiltonian and its fixed point behavior, 
the Luttinger li quid , . In the first subsection, we summarize well known results for the renormalization group flow of 



the parameters 



In the second subsection, we show how to apply perturbative techniques to the small cut-off 



theory to compute the uniform magnetic spin susceptibility. 



A. Renormalization group 

Following the Kadanoff- Wilson renormalization group procedure of Ref.0 we derive the one-dimensional scaling 
results that are essential for the calculation of the magnetic susceptibility. One starts with the partition function of 
the one-dimensional electron gas expressed in terms of a functional integral over anticommuting fields, namely 

Z= [ [ D^*D^e s[r ^ (9) 



The Euclidean action S = S° + Si corresponding to the g-ology Hamiltonian (Q), consists in a sum of free and 
interacting parts. Using the definition 

^ P Ak, ikn) = J% E f dre-^+^^Ar^r) (10) 
V L , =1 Jo 

with the fermionic Matsubara frequencies k n — (2n+l)7rT, n = 0, ±1,±2, ... and the notations k — [k,k n = 
(2n + 1)ttT] and q = (q, iq m = 2mirT) the two parts of the action are given by 
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and 



= E G^ 1 (fc)^, CT (fc)V P ,.(fc) (11) 



rp 



{ P,fel,fe2,9,o-} 



( fc i) ^*- P ^{ k 2)^~p^ 3 {k2 + q)ip P ,aA k i - v) 



2 1 _ 

{ p,fci,fe 2 ,g,(T,o''} 



To rewrite this equation, the gi term in Eq.(Q) has been subjected to the change of variable, q — > fe — fc' — 2pfci? — In 
the usual g-ology terminology described above, the constants 51 and g 2 stand for the backward and forward coupling 
constants between right and left moving carriers whereas 54 is the forward scattering amplitude for electrons on the 
same branch. The sums over wave vectors in Si[ip*,ip] are restricted by the bandwidth cut-off E = 2wi?fc - The 
measure is given by 

DrDiP = J] #; )CT (fe)# P)CT (fc) (13) 
while the one-dimensional free propagator of the fermion field has the form 

^ = -(^^M^) =ik^^ (14) 

with e p {k) as in Eq.Q. 

In the bandwidth cut-off scheme, the renormalization group procedure in one dimension consists in successive 
partial integrations of fermion degrees freedom in the outer band-momentum shell corresponding to energies ^Eq(1) > 
e p (k) > \Eq{1 + dl) for electrons and — \Eq{1) < e p (k) < — ^E (l + dl) for holes, with E (l) = E e~ l an effective 
bandwidth cut-off at step I > (Eq = 2Ep). For each partial summation, a complete Matsubara frequency sum is 
performed. Making use of the linked cluster theorem, the partial trace can be formally written as 



Here the tp s refer to fermion fields with a band wave vector in the outer momentum shell while the tp's pertain to 
lower momentum degrees of freedom that are kept fixed. The outer shell averages (connected diagrams) are defined 
by 

((T....?))o, c = Z- 1 j f rhpDQty...^) e s °^^ (16) 



From Eq. (|l5|) , the two- and four-point vertices of the Euclidean action are successively renormalized as a function of 
I. We can neglect the effect of 94 at this stage and reintroduce it later because it leads to logarithmic terms only 
in higher order, as explained at the end of this subsection. Taking into account the possible differences between 
parallel and anti-parallel spin scatterings, three distinct four-point vertices can renormalize. Upon renormalization 
one cannot tell the difference between gu and g 2 \\ - One is free to include all the corresponding scattering processes in 
a quantity defined as g±\\ (/ + dl) — g> 2 || [I + dl). Adding and subtracting then a coupling g 2 u (I + dl) = g 2 ± (I + dl) to 
the renormalized action, one is left with a form that under renormalization will preserve the structure of the original 
action, namely g±\\ (I + dl) — gi± (I + dl) and g 2 \\ (I + dl) = g 2 ± (I + dl) . In the notation of Ref.EU, one can write 



G 



/. 8 j"> i s, 

gjd+dl) = g(l) + ( 1 + £/ 



FIG. 2. Diagrams for the renormalization group calculation of the recursion relations for g\ and Q2- Wiggly lines are for 
interactions, dashed (solid) lines are for fermions on the left — ItF (right kp) branch of the spectrum. The bar on a line indicates 
that the corresponding integration is in the small momentum shell being integrated. 

z^(l)G° p (k) -> z^(l)z^(dl)G° p {k) (17) 

*2,3(J)ri,2 -» «2,3(0«2,3(di)ri,2 (18) 

where the renormalization factors Zi f2j a(dl) pick up the logarithmic outer shell contributions (0(dl)) of the cumulant 



expansion in Eq.(|l5|). If one performs the field rescaling z*ip(*' — > in 5 so that the free part of the Hamiltonian 
is invariant, one gets the recursion relations for the coupling constants 

gi,2(l + dl) = z 2<a {dl)zi 2 {dl) 9l , 2 {l) (19) 

At the one- loop level, there is no logarithmic one-particle self-energy corrections so that z\ = 1. The Hartree-Fock 
diagrams redefine the chemical potential in such a way that Luttinger's Fermi surfacecJ is preserved. As for the four- 
point vertex renormalization factors, z 2 (dl) and Zs(dl), they are obtained from the evaluation of the n — 2 outer-shell 
averages where 

S/f^*, ip,i/}*,ip] ~ ip* p ip*_ p ^-pip p + ip* p if}_ p ^*_ p 'iljp + ... (20) 

corresponding to outer-shell decompositions in the Cooper, and Peierls channels. An essential characteristic of the 
one-dimensional electron gas that emerges at the one- loop level is the quantum interference between both channels. 
Thus when the logarithmic diagrams are evaluated at zero external Peierls and Cooper variables, several diagrams 
cancel and only a 2Uf electron- hole bubble remains for the renormalization of g\ whereas g 2 is affected by a Cooper 
ladder graph. The relevant diagrams appear in Fig. 2. After the outer shell integration one is then left with, 

z 2 {dl) = l-g 1 {l)dl (21) 

i — 2 

z 3 (dl) = 1 - -^-dl, 

232 



where gi. 2 = gi.2/ (ttvf)- From the multiplicative renormalization (|19[), this gives the one-loop scaling equations 

d{2g 2 - gi) = 
dl 
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where the combination of coupling constants 2^2 — 5i is renormalization group invariant which is related to the 
conservation of the particles on each branch for g\ and g 2 scattering processes, when umklapp processes are neglected. 

Another property of interest is that the decoupling of the above two scaling equations can be understood as a 
consequence of the separation between spin and charge long-wavelength degrees of freedom. This is clearly manifest 
when the Hamiltonian representation Hi (I) of the effective interacting part of the action at step I is written in the 
Landau channel in a rotationaly invariant form 

H I [g 1 (l),g 2 (l)} = 5> 32 - gi ) Pp (q)p- P (-q) - 5i (0 £ S p (q) ■ S_ p (-g). (23) 

Here the operators for the long-wavelength charge (p p ) and spin (S p ) degrees degrees of freedom of branch p are 
respectively given, in operator form, by 

1 v — f 

p± P (±q) = — t=2_^ c ±p,*( k ) c ±pA k ± q) 



1 v-^v * f 

S± P (±g) = — -JjJ^ c ±pA k )va,f3C ±p j 3 (k±q). (24) 

with the vector a whose components are the usual Pauli matrices. Because of the integrated degrees of freedom, 
the summations on band wave vectors k are restricted to the interval | e p (k) |< Eq(1)/2. This way of writing the 
effective Hamiltonian together with the scaling equations (|2^) emphasizes that interfering and logarithmically singular 
correlations of the Peierls and Cooper channels do influence uniform correlations of the Landau channel. This must 
be taken into account in the calculation of the uniform magnetic susceptibility as we show shortly. 

Before concluding this subsection, we mention how 54 appears in the partial trace operation Eq. ([l5]) . One can show 
that the coupling of g± and 32 to 34 at higher order, also leads to logarithmic corrections for z 2 and 23, which are 
equivalent to a renormalization of the Fermi velocity in the scaling equations (p^). One then has g\ — gi/(jrv a ) and 
2<?2 - 91 = (2.92 - gi)/(nv p ), where 

v p , a =v F [l±g 4 (2irv F )- 1 }. (25) 
Here v a and v p are respectively the velocity of spin and charge excitations. Thus for g\ for example, one can write 

l + g 1 (Trv a ) H 

while the chatge coupling 2^2 — 9\ remains an invariant. As for the interaction term 34 itself, it does not renormalize 
at this ordenl3. 



B. Spin and charge susceptibilities from auxiliary fields. 

The renormalization described above is valid as long as the cut-off energy Eq(1)/2 is larger than the temperature. For 
smaller cut-off, the contributions from loop integrations are not logarithmic because of the Fermi occupation factors. 
While one could in principle modify the recursion relations to account for this and continue their integration^! it 
is simpler to apply a sharp cut-off procedure. We simply use the effective action obtained for I \sx{Ep/T) and 
apply perturbation theory to obtain the uniform susceptibility. We go through the exercise of formally generating 
the perturbative result through the auxiliary field (Hubbard-Stratonovich) method. This will ftUow us to exhibit 
the accuracy of standard perturbative techniques for a linearized spectrum in the Landau channeD'til and to find the 
infinite slope of the susceptibility in the zero temperature limit. But first, we must analyze the part of the Hamiltonian 
that we will take as the unperturbed one. 



1. The g4 theory 

Rescaling wave vectors to recover the original units, all band momenta |fc| are now smaller than fc ~ T/vp and, 
correspondingly, transfer momenta \q\ have to be of the same order. Hence we are now working with a small cut-off 
theory involving only low-frequency, low-momentum interactions within a single branch, 
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k a +pk F i 

1# = 4,* (*) + || E EE 4,* (*) 4,^ (*')<*,*'(*' + q)c p , a {k-q). (27) 

fc >°" ~ k,k' = -k +pk F 1 v,a' = \ 

It will be useful then, later, to work in an interaction representation where in the zeroth-order Hamiltonian the two 
branches ±fc.p do not interact with each other. We will need to know the value of averages such as (S pa (q) S p , a (<f )) 4 
computed with the above Hamiltonian. The label a refers to spatial direction, x,y,z, of the spin operator. Because 
the Hamiltonians H^ are small cut-off Hamiltonians, the exact irreducible vertex in the particle-hole channel can be 
taken as g± without fucjther Landau-theory-like renormalization from high-energy processes. A number of cancellations 
occur for this model,eZl so that the final expression has an RPA like form, 

I X° (q) 

(Spa (5) Sp> a (q )) 4 = S P ,p'S~_~,-- ! P ■ (28) 

I I — 294X P \q> 

In this expression, the non-interacting susceptibility x p (q) on one branch branch p = ± is given by 

x ; (a = _ 2E iMfiWk|±2)l, (29) 

p iq n + e p (k) - e p (k + q) 

where / (e p (fc)) is the Fermi function. At iq n = 0, this is one half of the total bare susceptibility since the sum over 
wave vectors is peaked near only one of the two Fermi points. At small wave vector and zero temperature, 

4(q,iq n ) = N(E F ) PVFQ . , (30) 
1 pvFq - iq n 

with N(Ep) = (ttvf)^ 1 the bare density of states per branch (half of the total bare density of states). Similarly, the 
charge-charge correlation function is given byEH 

i x° (q) 

(p P (q)Pp> (§0) 4 = — — 57^- ( 31 ) 
Note that in the iq n — 0, q — ► limit we have at low temperature, 

lun{S pa (q,0)S pa (-q,0)) 4 = - = -— (32) 

q^O 4 1 — (74/ (27TUf) 4 

km ( Pp (q, 0) Pp (q, 0)) 4 = 7 , /!T,l"l^ ^7^' ( 33 ) 



^O^vn, ,^vi: /,4 4 1 + 54 / (2ttw F ) 4 7TW 



where the spin and charge velocities are defined as above in Eq.(25) 



2. Auxiliary-field representation 

The magnetic susceptibility is obtained from a derivative with respect to an external magnetic field h. We choose 
units where the g factor times the Bohr magneton equals unity. In the presence of h, the partition function takes the 
form, 



Z{h}=Tr\e-e( H + +H -)T T exp 



- J^H I [g 1 (l),g 2 (l),T]-(^2^ P (Q,r) ) j h(-g,r)| dr 



(34) 



where we use the interaction representation in which the Hamiltonian Hp studied in the previous section plays the 
role of the unperturbed Hamiltonian. Using Gaussian integration (Hubbard-Stratonovich decomposition) , to decouple 
the interactions between branches, this may be rewritten as 

Z = ZiJ J V(j)VM. exp ^- EE dT IMQi T^-pi-q, r) + M,(j, r) ■ M_,(- g , r)]j 
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x ( T T cxp - / drH [</>, M, h, t) 



(35) 



where <fi an d M are real auxiliary fields for charge and spin degrees of freedom respectively and 

W[0,M,h,T] = E [ 2 W 2 92 - 91 P P (q,r)4> p (-q,T) 



+ 2V9i(0 S p (q, t) ■ M p (-q, r) - S p (g, r) ■ h(-q, r) 



(36) 



while Z4 — Tr e > 3 ( H ++ H -) an d (...) 4 are respectively the partition function and corresponding averages calculated 

with in the presence of the 54 interaction only. As in previous sections, we use the following definition for quantities 
in Matsubara frequencies uj rn — 27rmT, m — 0, ±1, ±2, ... 



h(r) = Vr£V*"» r h(w ro ) 

h (u> m ) = VT / dre^ mT h (r) . 
Jo 



(37) 
(38) 



With these definitions, and the change of variable M p (q) — > M p (q) + h (q) / , q = (<?, w m = 27rmT) the 

zero-field dynamic magnetic susceptibility per unit length can be formally expressed in terms of an average over the 
magnetic auxiliary field M, that is 



S 2 lnZ[h] 



5h a (-q)Sh a (q) 



h=0 



5i (0 



(39) 



where the subscript a stands for the orientation of the magnetic field. 

Applying the linked-cluster theorem in Eq.([35|) and evaluating the averages 
previous subsection allows one to write, in the rotationally invariant case, 



using Eqs.(B3) and (p3|) of the 



Z = Zij J V</>VMexp J - E [<Pp{Q) A P,AQ)M-q) +M p (q)B p , pl (q)M p ,(-q)] +O (M 4 ,0 4 ,MV) 



with the matrix elements 



a p,p&) = 2(292 ~ 91) 



X° P (0) 



2 94X P (?) 2 



(252 - 9i)x°c,p (5) 



(40) 



(41) 



A -,+ = 1 



and 



BpJq) = --9x (I) 



xl (5) 



2 JiW 1-|94X°(§) 



9i (0x°, P (§) 



(42) 



(43) 



= = 1 



(44) 



respectively for the charge and spin degrees of freedom of the Landau channel. 

In these expressions, g\ (I) is given by Eq.(p6|) with I sa ln(A/T) where A is a cut-off of the order of the Fermi 
energy, hence, 
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The neglect of mode-mode coupling or anharmonic terms in Eq. (J4C|) is quite justified in one-dimension. Indeed we 
are in a low cut-off theory so that only the q = components of the Hubbard-Stratonovich fields are important. When 
the unperturbed Hamiltonian is quadratic in fermion fields, the coefficients of the q — mode-mode coupling terms 
vanish because they are derivatives of the density of statespi a quantity that is a constant for a linear dispersion 
relation. For the same reason, the spatial rigidity of correlations of the auxiliary fields cf> and M vanishes. Here the 
unperturbed part of the Hamiltonian is the interacting g 4 theory. It is known that in the thooty with a linearized 
dispersion relation, the density fluctuations are Gaussian, hence there is no mode-coupling term.cll Furthermore, there 
is no singularity in this theoryO. Here we will take into account the fact that the dispersion relation is not linear. 
This could make mode-mode coupling terms become different from zero at high temperature, far from the Luttinger 
liquid fixed pointtS. 



3. Susceptibilities 



The Gaussian fluctuations of M evaluated with the functional Eq. 
Eq.(B9f). We find in the rotationaly invariant case 



give us the magnetic susceptibility through 



x(i) 



(46) 



Using the expression Eq. ([43|) for B p p , one obtains 



x(q) 



1 

4 



91 (0 X°,+ (q) X° (0) + (x° + (?) + 5&- (§)) 



HI (0x2,+ (M,-®-i 



In the Luttinger-liquid limit, (linear dispersion relation) one finds for the retarded spin susceptibility, 



1 2 



(v s qf 



4 irv s (a, + irjY - ( v s v s ) q 2 



u = v a 



v =v rr - 



gi(0 

9i{l) 



vf — -z— 



.94 , gi (I) 



2tt 2tt 



94 gi(l) 



(47) 

(48) 

(49) 
(50) 



2tt 2tt 2tt 

which reduces to the known result@,E3 in the limit g± = 0. 

For the static susceptibility of interest to us, the result is simpler. In that case, we take u> — > first and notice that 



Xl+ (q = 0,uj = 0) = xl,- (q = 1 uJ = 0) = X ° (T) , 
where gi (I) is evaluated at Ij- = h\(Ep /T). This leads to 



X(X) 



1 ( X° P (T) 

2 \l-±g4 X « 



(T) 



(51) 



(52) 



To leading order in g\ (T) , one finds 



with 



x(T)«^(T)[l + iff 1 (T) X °(T) + „.] 



X° P (T) 



1 - (2ttw f ) 1 g 4 



(53) 
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In the absence of </4, this expression, coincides with the result obtained previously by Dzyaloshinskii and Larkinta and 
by Lee et alB. Our more general expression for the susceptibility Eq.(|52]) may be rewritten as 

x{T) = — w+nvwm' (54) 

For the repulsive sectcp-|<?i > 0, gi(T — > 0) — * is irrelevant in the zero temperature limit and there one recovers the 
Luttinger-liquid resultEi, 

X (T = 0) = ( £ S pa (0, 0) 5 p , a (0, 0) ) = — r =I - 1 =-^— (55) 

I 2tTVf[1 — (ZTTVf) 34j 2TTV a 

The analogous calculations for the charge fluctuations also lead to similar results. It suffices to do the following 
substitutions in any of the above spin susceptibility results: 

.94 -> -54 (56) 

.91 ->■ gi - 2g 2 (57) 

In particular, in the static limit, the charge susceptibility Xc (T) (or equivalently the isothermal compressibility 
K T (T) = Xc{T)/n 2 ) is given by 

Xc (T) = ( £ Pp (0, 0) Pp (0, 0) \ = — - J X ° P — — . (58) 

\rri / 1 + ^ .94 + 2o 2 - gi )r„ (T) 



\p,p' 



Again in the repulsive sector g\ > 0, 2g 2 — g\ is a renormalization group invariant when umklapp-scattering can be 
neglected (cf. Eq. and one recovers the Luttinger-liquid result in the zero-temperature limitE3, 

Xc {T = 0)= 2Trv F [l + (27rv F )-i(g 4 + 2g 2 - gi )} (59) 



4- Infinite slope in the zero temperature limit and third law of thermodynamics 
While in the absence of umklapp scattering the charge susceptibility comes in the zero-temperature limit with zero 



slope, as in a Fermi liquid, the dependence of the magnetic susceptibility Eq.(54) on the marginally irrelevant variable 
gi (T) implies an infinite slope in the zero temperature limit. This can be seen as follows. Since dXp/dT has a finite 
limit as T — > 0, the singular contributions comes from 

9X(T) ^ q (t MT) 
dT x 1 ' dT 

The temperature derivative of the marginally irrelevant variable equals infinity at T = as can easily be obtained 
from the temperature derivative of Eq. (|45| ) 

dgi(T) _ £1 £i_ 1_ 

hence, 

lim - = 00. (61) 

A superficial look at the zero temperature infinite slope Eq.(|6l|) suggests a violation of the third law of ther- 
modynamics. Indeed, consider the grand potential Q = E — TS — fj,N, whose differential change is given by 
dfl = —SdT — Ndfi — M-dB. Normalizing to unit volume on finds, 
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T=0 



d 8M 
&T~dB 
d 2 



dB 2 dT 



n 



d d 2 VL 
dfdB 2 
d 2 S 



T=0 



T=0 



dB 2 



(62) 
(63) 



T=0 



The infinite second derivatives of the entropy seem to contradict the third law of thermodynamics that says that 
the entropy at zero temperature is independent of external parameters. However, because the point T — 0, B = is a 
critical point (infinite correlation length), the free energy at this point is not analytic and we are not allowed to invert 
the order xtf differentiation as we did on the second line. Hence, one cannot conclude that ^x| T _ = 00 v i°l^t es the 
third lawBJ The entropy of the Luttinger liquid does vanish at zero temperature, independently of B and fi. 



V. COMPARISONS WITH MONTE CARLO SIMULATIONS 



In this section, we shall compare the RG results first with the zero-temperature exact results and then with all 
the finite-temperature results exhibited in Fig.l. It is not yet possible to do simulations at low enough temperature 
to confirm or not the existence of the infinite slope seen in the RG approach in the T — > limit since this also 
requires huge system sizes. It should however be possible to verify that regular extrapolation of finite T results to the 
T = limit is not possible. At higher temperature and larger interaction strengths, non-logarithmic terms become 
important. We will show that it is possible to estimate these. In two dimensions, it has been possible to explain-tha 
complete temperature-dependent magnetic properties of the Hubbar.d model by using diagrammatic approachest2l,E2l 
far from half-filling or the Two-Particle Self-Consistent approacho,t3 at arbitrary fillings. We briefly explore the 
predictions of these approaches in one-dimension. The reason for their failure in one dimension will help understand 
the correct way to proceed. 



A. Comparisons with the Renormalization Group approach 

To compare numerical Monte Carlo results with our RG results for the susceptibility Eq. ( p2[ ) or equivalently Eq. ( p4| ) , 
(remember the trivial factor of four to compare with simulations) one needs to know the initial values of g\ and g± 
entering the scaling equations. As argued above, it is not strictly correct for the Hubbard model to assume that these 
constants can be taken as gi = g± = U because the dispersion relation is not linear, and the cut-off in the initial model 
is not as in the g-ology model. Let us start by a comparison with the zero temperature exact results of Shiba. We 
find that our result Eq.(|5^) differs from that of Shiba by at most 15% up to U = 4 if we choosey = U. At U = 2, the 
RG result is accurate to about 1%. The same conclusions were reached a long time ago in Ref.E3. Hence, we conclude 
that at small coupling the estimate g\ = 94 = U should be accurate. 

To do a more general comparison at finite temperature, our results Eq.(|54|) and ( [i"5"|) require in addition a value 
for the cut-off A. Taking A = 2, which corresponds roughly to half the bandwidth, of the order of the Fermi energy, 
produces the results in Fig. 3. At low temperature, the results are not too sensitive to the value of A which enters 
only logarithmically. The comparison in Fig. 3 shows that the disagreement with Monte Carlo results does become 
larger as U increases. Higher order logarithmic terms will not change the picture since they are smaller. As shown in 
the following subsections, the main source of discrepancy resides in non-logarithmic contributions that are actually 
not negligible at high temperature where the non-divergent (non-logarithmic) contributions dominate the magnetic 
susceptibility. 

Le us discuss in turn the various features appearing in Fig. 3. The Monte Carlo data shows that the temperature at 
which the maximum susceptibility occurs is essentially independent of the interactions. The interactions determine 
only the overall enhancement and sharpness of the maximum but not the position. The position of the maximum 
T max within RG also does not depend strongly on interactions but there can be a 10% shift in T max when one changes 
the cut-off from A = 1 to A = 2. However, since the RG results do not seem very reliable near the maximum, we 
postpone this discussion to Sec.V.E. 

Decreasing the temperature from the position of the maximum, the RG predicts a first inflection point located at 
Ti = 0.099 ±0.002. This temperature is just at the limit of our Monte Carlo data. The inflection point can also not be 
seen in the more recent data of Ref.cl since only the points T = 0.1 and T = 0.05 are available in the low-temperature 
region (See Fig.l). This inflection point is again an effect that is caused by the band structure and within the above 
accuracy the position is independent of interactions. Indeed, while the curvature of the non-interacting susceptibility 
d 2 xo/dT 2 is negative at the maximum, it must be positive at low temperature because of the characteristics of the 
one-dimensional band with a parabolic bottom. The position of this inflection point is the same, within the quoted 
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Fig. 3 




0.00 0.10 0.20 0.30 0.40 

Temperature in units of t 

FIG. 3. Magnetic susceptibility \ evaluated from Eqs.^5^) or equivalently Eq.j54|) with the naive replacements g<t = 
gi — U and, using Eq.(Q), gi(T) = g\j (l + In ~J with A = 2. The spin velocity is obtained from Eq,(pi|), 
Vp,<y = Vf[1 ± <?4(27t?;f) _1 ]. The inset shows the low temperature region. In the inset, the X symbols show the low-temperature 
limit predicted by the RG. 

accuracy, whether one uses A = lorA = 2ora temperature independent interaction as in Sec.V.E. With a parabolic 
band instead of a cosine band, one finds an inflection point at Tj = 0.101 ± 0.002, again essentially independent of 
cut-off or interactions. 

Since the zero temperature value of the susceptibility is reached from finite temperature with an infinite slope, 
as discussed in the previous section, there is a second low-temperature inflection point. For A = 2, it appears at 
Tl = 0.0329, 0.0328, 0.0312 for, respectively, U — 2,3,4. This temperature is roughly three times lower than our 
lowest temperature in the Monte Carlo simulations and also lower than the other lowest available numerical results 
exhibited in Fig.l. The infinite slope in the susceptibility predicted by the RG appears confined to a small temperature 
range even when seen on a magnified scale, as one can check in the inset of Fig. 3. 

Contrary to the previous T max and Tj, identified above, the location of the low-temperature inflection point Tl is 
not purely a band structure effect. It is clearly absent from xo (T) ■ The low-temperature inflection point appears in 
the RG calculation because of the competition between the band-structure effects that lead to d 2 \o/9T 2 > and 
the logarithmic singularity in g\ (T) that leads to negative curvature in the full susceptibility at low temperature. 
As long as the interaction is sufficiently strong, this low-temperature inflection point occurs at a temperature that is 
remarkably independent of interactions. For lower values of the interaction, U — 0.5 and {7 = 1, one finds that the 
inflection point is at Tl — 0.0208 and Tl = 0.0288. Comparing with the above results, one sees that for 1 < U < 4, 
the low-temperature inflection point is independent of interaction within about 10% (Tl = 0.031 ± 0.002) while 
for 0.5 < U < 4 the position varies from 0.02 to 0.03. The weak dependence on interaction comes from several 
factors: a) The interactions influence the susceptibility weakly, mainly to first order in the expansion of the RPA-like 
denominators, b) The inflection point occurs in a regime where g\ (T) takes its asymptotic form, ~Kv a f In (A/T) which 
depends on the interactions only weakly through v a . c) Only the order T 2 of the Sommerfeld expansion of \o suffices 
to obtain an accurate result. These three simplifications allow one to obtain an accurate analytical expression for the 
location of the inflection point. However, since terms up to order T 2 ln 3 (A/T) must be kept, the resulting equation 
is transcendental and must be solved numerically. 

The location of the low-temperature inflection point Tl clearly depends on the value of the cut-off A. However, 
since the dependence is logarithmic, the above results are not so sensitive to the precise value of the cut-off. For 
example, for A = 1 one finds T L = 0.0219, 0.0304, 0.0349, 0.0345, 0.0324 for, respectively, U = 0.5,1,2,3,4. These 
values of Tl are larger than the corresponding values for A = 2 by at most 6%. For a parabolic band and A = 2, one 
finds Tl = 0.035, 0.035, 0.034 for U — 2,3,4 values that are at most 8% larger than the corresponding results for the 
cosine band with the same A = 2. 
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B. Diagrammatic Kanamori approach 



From now on, we concentrate on the susceptibility at temperatures above the low-temperature inflection point, 
where the RG singularities do not contribute appreciably to the susceptibility. It was suggested long ago by Kapamori 
that the interaction appearing in RPA expressions should be renormalized. Following this idea, it was shownlij that 
in two dimensions the renormalized U rn can be accurately computed as follows: 

Urn = ( i+t^ib^o) ) ■ (64) 

In this expression, the quantity Xpp (QMn = 0) is the Cooper bubble for a total incident momentum Q. The average 
is over the values of Q. Several different types of averages can be done, namely over the whole Brillouin zonejla or 
over wave vectors corresponding to Wues of Q = k + k' such that both k and k' are within an energy equal to tha 
temperature T of the Fermi surface.L3 The latter type of approximation is closer in spirit to Fermi-liquid theoryEZl 
and gives overall better results. 

When the one-dimensional version of the procedure of Ref.tS is applied to the present one-dimensional case, the 
agreement with the Monte Carlo results is apparently much better at high temperature than with the scaling approach. 
The temperature dependence is essentially correct and there is simply an underestimation of the overall magnitude 
of the susceptibility by 3%, 4%, and 9% for, respectively, U = 2,3, and 4. The value of U rn that we find may be 
approximated by 

U rn « ^ (65) 

with A « ( Xpp (QMn = 0)) « 0.29. 

In two dimensions, the agreement with Monte Carlo simulations is much better than found here and it is valid for 
all wave vectors, as long as the filling is such that there is no zero-temperature phase transition. In one dimension, 
the bare susceptibility at q = 2k f diverges logarithmically with temperature at any filling because of nesting. This 
means that at sufficiently low temperature, RPA with a temperature and wave vector independent U rn predicts a 
finite temperature phase transition at any filling. This is prohibited by the Mermin- Wagner theorem. Sure enough, 
for the values of U studied here, this transition occurs at a temperature much lower than those investigated above 
with Monte Carlo, but nevertheless, this is a question of principle that cannot be overlooked. 



C. Two-particle self-consistent approach (TPSC) 

The TPSC approacbH-H3 avoids any finite-temperature phase transition in both one and two-dimensions. Hence, we 
may check wether this gives a better agreement with Monte Carlo simulations than the previous approach. The one- 
dimensional version of the theory can be summarized as follows. One approximates spin and charge susceptibilities 
Xspj Xch by RPA- like forms but with two different effective interactions U sp and U c h which are then determined 
self-consistently using sum rules. Although the susceptibilities have an RPA functional form, the physical properties 
of the theory are very different from RPA because of the self-consistency conditions on U sp and U c h- The necessity to 
have two different effective interactions for spin and for charge is dictated by the Pauli exclusion principle (n^) — {n a ) 
which implies that both Xsp and Xch are related to only one local pair correlation function {n^n^). Indeed, using the 
fluctuation-dissipation theorem in Matsubara formalism and the Pauli principle one can write: 



l&„iS=,,-2(„» l ) = J_^_|gL_, (67) 

9 9 

where (3 = 1/T, n = (ri|) + (ni), q= (q,iq n ) with q the wave vectors of an N site lattice, iq n Matsubara frequencies 
and Xo(q) the susceptibility for non- interacting electrons. The first equalities in each of the above equations is an exact 
sum-rule, while the last equalities define the TP-SC approximation for Xch{q) and for Xsp(o)- ln this approach, the 
value of (rif ^x) ma y be obtained self-consistentlya by adding to the above set of equations the relation U sp = 3|x(0) U 
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with <?tj.(0) = (n|n|)/(n|)(ri|). As shown in Ref.E-3, the above procedure reproduces both the Kanamori-Brueckner 
screening described in the previous section as well as the effect of Mermin- Wagner thermal fluctuations, giving a phase 
transition only at zero-temperature in two dimensions. In two dimensions, there is however a crossover temperature 
Tx below which the magnetic correlation length £ can grow exponentially. Quantitative agreement with Monte Carlo 
simulations is obtainedo for all fillings and temperatures in the weak to intermediate coupling regime U < 8t. The 
equation for charge Eq.(^) is not necessary to obtain the spin structure factor. 

In one dimension, the absence of phase transition in this theory at finite temperature can be proven as follows. Near 
the temperature at which the phase transition would occur in RPA, SU = U m f, c — U sp ~ 0, (U m f. c = 2/xo (2kp,0)) 
the q = 2kp susceptibility at zero Matsubara frequency is becoming very large so that the self-consistency relation 
Eq. (p7|) can be approximated by 

2r /S t /, A ^ + ^ ) -"- 2 <""">- c - < 68 > 



where 



f 2 = 
?o — 



1 d 2 X o(q,0) 



2xo (<?) d 2 q 



(69) 

q=2k F 



^^(U^/SU) 1 / 2 , (70) 

and where the integral is for q around 2kp or —2kp and C contains contributions from non-zero Matsubara frequencies 
and from corrections to the Lorentzian approximation used for the iq n = contribution. Then, 

iT [<hl ' - 4T ^ = n-2(n 1 n l )+C. (71) 



U sp ej 2n^- 2 +q 2 U SP Z 2 

Since the right-hand side is a finite number, £ behaves as U sp ^ 2 / (4T), becoming infinite at zero temperature only. 
There is however nothing to prevent a zero-temperature phase transition in the theory, so that it cannot describe 
accurately the one-dimensional systems at very low temperature. 

Comparisons with the Monte Carlo simulations reveal discrepancies of order 25% for the case U = 4. This means that 
contrary to the two-dimensional case, this approach does not even reproduce the Kanamori-Brueckner result described 
in the previous section. If we had taken (riffii) from Monte Carlo data in the sum rule Eq.(|67|) instead of computing 
it self-consistently from the ansatz U sp — U {n^n^) / {n^){rL^) the results would have been much better. In other 
words, the calculation of (nf7i|) self-consistently is an assumption that fails in one dimension even more drastically 
than the RPA functional form with an effective U sp . The whole approach fails completely at low temperature. 
Indeed U m f. c = 2/\o (2kp,0) vanishes as T — > while 5U = U m f, c — U sp has to remain positive according to the 
above arguments. This in turn implies that U sp tends to zero at zero temperature which means that the uniform 
magnetic susceptibility would not be enhanced at zero-temperature with this theory, contrary to both exact and 
renormalization group results. One could have hoped to use this approach to evaluate non-logarithmic contributions 
at high temperature and inject them in the RG expression, but in fact the logarithmic temperature dependence of 
the susceptibilities starts at rather high temperature. 



D. General reason for the failure of higher-dimensional approaches 

All of the above approaches fail for several reasons. They do not take into account the destructive interference 
between Cooper and Peicrls channel and they try to describe the whole q dependence of the susceptibility with a 
single wave- vector-independent effective interaction U sp or U rn . This is very different from the scaling theory which 
clearly shows that the effective interaction near q = is different from that near q = 2kp. 

The fact that the magnetic susceptibility in one-dimension cannot be described with a single w*we vector independent 
effective interaction U sp or U rn is illustrated in Fig. 4a. The simulations are taken from Ref.(l23). The solid lines are 
for an RPA-like form 

xoM) (72) 



and the data is normalized by XofeO) as in Ref.(E3). While the bare value for the simulation is U = 2, one needs 
a renormalized value U = 1.3 to fit the components near q = 0, while to fit near q = 2kp one needs U = 1.7. Note 
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FIG. 4. a) Magnetic susceptibility \ (q) as a function of wave vector, normalized to its non-interacting q — value. The 
Monte Carlo data is from Ref.(E3). The solid lines illustrate that when the whole q dependence is considered, RPA fits are 
inadequate, even with renormalized values p£ the interaction, b) Magnetic structure factor Eq.(^) as a function of wave vector, 
S (q) . The Monte Carlo data is from Ref.(c3). The solid line illustrates that this time, as opposed to part (a), a simple RPA 
fit with a renormalized value of U may work misleadingly well. 

also that the simulations are done for the canonical ensemble so that the q — component is strictly zero and is not 
shown. 

The magnetic structure factor contains not only the above susceptibility, but also all the non-zero Matsubara 
frequency components of the susceptibility. That iSpwhy it is less sensitive to the ±2kp effects of one dimension. As 
illustrated in Fig. 4b, the Monte Carlo data of Ref.(ti3) for U = 2 can this time be fitted, misleadingly, with a single 
renormalized U = 1.5. 



E. Modified-Kanamori approach in one dimension 

In the renormalization group description of magnetic fluctuations, it is clear that the effective interactions near 
q = and near q = 2kp are different. No finite or zero-temperature phase transition occurs even though the uniform 
magnetic susceptibility is enhanced. However, this approach takes into account only logarithmic terms and applies 
only in the vicinity of cither q = or q = 2kp. It does not allow one to draw the full g-dependent susceptibility 
appearing in Figs. 4 for example. Furthermore, we have seen that the renormalization group result for q = becomes 
inaccurate compared with Monte Carlo simulations at high temperature. This can be understood as follows. In 
addition to the fact that we have employed an approximate sharp cut-off procedure, we do not know the exact value 
of <74 and of g± entering the recursion relations. But more importantly, logarithmic terms are less singular at high 
temperature so that non-logarithmic terms also become important. To estimate non-logarithmic contributions, we 
proceed as follows. 

In the same spirit as the Kanamori approach in higher dimensions, we want to find an effective interaction that 
contains the effect of other channels. This time however, the effective interaction should be valid only for q ~ 0. The 
RG result suggests that it is the contribution from the particle-hole channel that is important. Indeed, g\ entering the 
) for the susceptibility contains the set of all diagrams generated by summing 2kp electron- hole 
bubbles.Eil This can be seen from Fig. 2 and from the recursion relation Eq.(|2^). In Fig. 2, all cross terms involving g\g<i 
cancel each other in the cut-off theory so that only the g\ contribution represented by the bubble is left. This shows 
that in the RG approach only 2k f electron- hole contributions are important. In the computation of the magnetic 
susceptibility however, one can observe that internal summations over all momentum transfers are present (including 
q = 0_where electron-hole bubbles are also the only type of diagrams that contribute in the calculation of x f° r Hp 
Rcfs.(li3'EZl)). In order to include these non-logarithmic effects in the calculation of the susceptibility we average the 
resummed series of electron-hole bubbles over the entire Brillouin zone. More specifically, we take 



expressions (|52|) 



and ( 54 
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FIG. 5. Numerical results for the temperature dependent susceptibility, as already illustrated in Fig.l, compared with the 
Modified-Kanamori approach defined by an RPA from with a renormalized value of the interaction U m given as a function of 
the bare value U by Eqs.Q or @. 

Um = (i + u Xeh U ( q ^ qn = o)) q > (73) 

with 



:h (q, iq n = 0) = —N 1 



. e {k-q)-e{k) 



Here the full effect of the lattice is restored by taking e (k) — 2t cos ka with the summation over k that covers the 
entire Brillouin zone [— — ,—]■ It is clear from this procedure that for \q\ « 0, the average will contribute to a g& type 
of process while for | q |w 2kp the contribution will be to a g\ type of process, as suggested by Eq.(^). One can 
also check that the Cooper and Peierls-type of bubble diagram, such as those appearing in 31172 processes in Fig. 2, 
do cancel at half-filling when averaged over the Brillouin zone. At this filling, there is perfect particle-hole symmetry, 
even with a 2t cos ka dispersion relation. Even if this cancellation is no longer strictly valid away from half-filling, we 
write 

X(T) = - - (74) 

for the total susceptibility with U m computed from Eq.(|73|) 

The results for the susceptibility are compared with the Monte Carlo results in Fig. 5. The agreement is within the 
statistical uncertainty except near T = 0.33. As in the Kanamori approach, for the range of temperatures illustrated 
in Fig. 5 we can approximate U m by a temperature-independent value given by 

u - " T^uk, < 75 > 

with A m = (A m (q, iq n — 0)) . This last approximation gives, within about 2% in the worse case, the same value for the 
susceptibility as that obtained from Eq.(f73|). Using the low temperature value, A m ~ 0.25, over the full temperature 
range improves the agreement near T = 0.33. At T — the above modified-Kanamori approach overestimates the 
exact result by 13% at U = 2 and by as much as 17% for U = 4 showing clearly that a Fermi-liquid like extrapolation 
of the finite-temperature data does not yield the correct zero-temperature limit, as should be clear from Fig. 5. 
Independently of the RG then, it seems_certain that a change in curvature is needed to extrapolate to the correct 
zero-temperature limit. The RG resullo with ^-.U on the other hand does have a change in curvature below 
T = 0.1 and extrapolates to the exact T — resulted to within 1% for U — 2 and makes physical sense even if 15% 
deviations occur in the extrapolation around [7 = 4. 
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Within the approach we just discussed, we can now come back to the question of the location of the maximum in 
the spin susceptibility. The Monte Carlo data suggests that the location of the maximum is essentially independent 
of interaction, a result that is an obvious consequence of the analytical form of the susceptibility that we used when 
U rn has a temperature dependence that can be neglected. We find that T max = 0.178 is determined purely by the 
band structure and at quarter filling is even the same at the 1% level for a cosine or a parabolic band. Nevertheless, 
we do not have a simple analytical expression for the position of the maximum. Indeed, a Sommerfcld expansion of 
Xo up to order T 6 predicts that xo always increases with temperature. It is the competition of this increase with the 
eventual Curie- Weiss decrease of the susceptibility at high temperature that produces the maximum. The Curie- Weiss 
decrease is beyond the radius of convergence of the Sommerfcld expansion. 

To extract the strength of interactions from a measurement of x (T) one would need to know the bare value of xo {T) ■ 
When xo (T) is unknown, one measure of the strength of interactions that suggests itself, given the position of the 
maximum T max and of the high-temperature inflection point Tj is R = (x (T max ) — \ (T)) Ix (Imax) • In percentage, 
one finds, for U = 2 and U — 4, respectively, R = 4.2% and 5.2%. Unfortunately it turns out that this result is too 
sensitive to the band structure to really be useful in practice since for a parabolic band the above cosine band results 
are replaced by R = 3.7% and 4.5%. 

VI. DISCUSSION AND CONCLUSION 

We have shown in this paper that the temperature-dependent magnetic spin susceptibility of the Hubbard model 
at quarter- filling has the following general features in the weak to intermediate coupling regime (0 < U < At). 
As temperature is decreased, one encounters a maximum in the susceptibility at T max ~ 0.18i that arises from 
the competition between the Curie- Weiss high-temperature decrease and the low temperature increase caused by 
the proximity to the characteristic van Hove singularity of one-dimensional systems with a parabolic band bottom. 
Decreasing the temperature, there is then an inflection point at <~ O.lt. The position of the maximum and of 
the inflection point are essentially independent of interaction strength, as confirmed by Monte Carlo simulations. 
Although a rough measure of the strength of the interactions may be obtained by comparing the relative size of the 
susceptibility at the maximum and at the inflection point this is not very reliable since one needs high accuracy as 
well as rather detailed information on the band structure: Indeed, a ratio (x (T max ) — \ (Tj)) j\ {T max ) of order 4% 
may correspond to either U ~ At or U ~ 2t depending on whether the band structure is parabolic or cosinusoidal. 
As shown in Fig. 5, the magnetic susceptibility curve obtained from Monte Carlo simulations between T ~ O.lt 
and T <~ 0.33i may be reproduced very simply by an RPA-like form with an effective interaction U m for the q = 
component of the susceptibility. For the case considered here, one finds U m = U/ (1 + U\ m ) with A TO ~ 0.25. That 
effective interaction may be computed from the bare one using the subset of diagrams that naturally appears in the 
RG calculation. This effective interaction allows one to take into account non-logarithmic corrections that are beyond 
the RG approach. The effective interaction is valid only near q = 0. Approaches to the interacting problem that work 
in dimensions larger than d — 1 fail essentially because they keep an RPA form for all values of q, an approximation 
that is incorrect in one-dimension because of the destructive interference that occurs between q = 2fcp Peierls and 
q = Cooper channels. 

At T = 0, the logarithmic terms in the leading marginally irrelevant operator of the RG lead to an infinite slope 
dx/dT. This Luttinger-liquid result is strikingly different from the Fermi liquid prediction d\/dT = at T = 0. Fig. 5 
clearly shows that a Fermi-liquid extrapolation of the finite-temperature data to the T = exact result is inappropriate. 
The RG approach that we used also allows to show that in the presence of band structure effects (parabolic band 
bottom) the appearance of the Luttinger-liquid regime at low temperature is signaled by a low-temperature inflection 
point located at Tl = (0.031 ± .002)i for t < U < At and A = 2t. Although the near independence of Tl on the precise 
value of the interactions for sufficiently large bare U does not have a simple origin, as discussed in Sec.V A, it comes 
in large part from the fact that in this case Tl occurs in a regime where the marginally irrelevant interaction takes its 
asymptotic form irv a / In (A/T) . This form depends on interactions only weakly through v a . The location of Tl is also 
not so sensitive to the (unkown) value of the cut-off since for A = t one finds that Tl is larger than the corresponding 
values for A = 2t by at most 6% over the whole range 0.5t < U < At. For A = 2t, Tl can also be larger by at most 8% 
for a parabolic band compared with a cosine band. The above Luttinger-liquid regime (where d 2 x/dT 2 < 0) occurs 
at a temperature lower than what has been achieved by numerical calculations up to now. It is quite a challenge 
to verify them. It is also important to realize that even though the logarithmic Luttinger-liquid limit shows up in 
the spin susceptibility at quite low temperature, logarithmic contributions do appear at higher temperature in other 
quantities such as the longitudinal spin relaxation time T\. 

Although our results are for the quarter-filled model, it is clear that all the qualitative features should hold for 
fillings that are not too close to half-filling. A good quantitative estimate for the location of the characteristic features 
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of the temperature-dependent spin susceptibility can easily be obtained for any other filling using the simple analytical 
expressions that we have found in Sees. IV. B and V.E. Our analytical expressions may also be used to compute the 
charge susceptibility (compressibility), which as we saw has no singularity in the T — > limit, as long as umklapp 
scattering can be neglectecS. One could also use our approach to make quantitative predictions for more general 
models than the Hubbard model. Comparisons with experiment should help to establish appropriate bare microscopic 
parameters for the Hamiltonians of one-dimensional organic conductors. 
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